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We study how macroscopic superpositions of coherent states produced by the nondissipative 
dynamics of binary mixtures of ultracold atoms are affected by atom losses. We identify different 
decoherence scenarios for symmetric or asymmetric loss rates and interaction energies in the two 
modes. In the symmetric case the quantum coherence in the superposition is lost after a single 
loss event. By tuning appropriately the energies we show that the superposition can be protected, 
leading to quantum correlations useful for atom interferometry even after many loss events. 
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In large systems, macroscopic superpositions of quan- 
tum states have extremely small decoherence times, mak- 
ing them impossible to observe [ij. The largest super- 
positions generated so far are superpositions of coher- 
ent states (CS) of a cavity field, and their progressive 
transformation into statistical mixtures as time evolves 
has been observed The first proposal to generate a 
macroscopic superposition of CS (MSGS) with light con- 
sisted in sending photons through a medium presenting 
a strong Kerr non-linearity In such media, the dy- 
namical phases of Fock states are nonlinear in the pho- 
ton number, thus the phase of an initial CS is split. In 
metastable vapors of ultracold bosonic atoms, interac- 
tions between atoms lead to similar nonlinearities, whose 
strength can be tuned experimentally by using Feshbach 
resonances 0, By trapping optically the condensed 
atoms in a double- well potential one realizes an external 
Bose-Josephson junction (BJJ); an internal BJJ is ob- 
tained by trapping in a single well atoms in two distinct 
hyperfine states coupled by a resonant field. In analogy 
with fight, the nonlinear dynamics generates MSGSs af- 
ter a sudden quench of the tunnel amplitude (for external 
BJJs) or a switch off of the coupling field (for internal 
BJJs) d, 01. Whereas MSGSs of fight are presumably 
destroyed after a single photon loss Q the situation is 
less clear for atoms. Up to now only squeezed states - 
which are produced at earlier times than MSGSs - have 
been observed in BJJs with a few hundreds of atoms 0- 
[lo| . In order to know if future experiments could produce 
MSGSs in BJJs, it is desirable to study how robust they 
are with respect to decoherence. 



Decoherence effects due to atom losses and phase noise 
on squeezed states have been analyzed in detail [ll- 13|. 
For what concerns MSGSs, only decoherence caused by 



noise due to photon scattering [l3| and phase noise (l5| 
has been studied so far. Under current experimental con- 
ditions the first noise is negligible and the second one can 
be reduced by using spin-echo (lol. Il6|: in contrast, atom 
losses are unavoidable. 

In this letter, we investigate how MSGSs generated in 
BJJs are affected by atom losses. Even though our anal- 
ysis also applies to one- and three-body losses [l3 we 
focus on two-body losses, due to scattering of two atoms 
in the magnetic trap which changes their spin and gives 
them enough kinetic energy to be ejected from the trap. 
The latter processes are particularly detrimental in the 
experiments of Refs. We first analyze the dynam- 

ics of a lossy BJJ from the point of view of state condi- 
tioning, assuming that the total number N of condensed 
atoms can be measured precisely, both initially and at the 
MSGS formation time tq. We study how much coherence 
is destroyed by a single loss event occurring at a random 
time between and tg. We find quite different answers 
depending on the degree of asymmetry between the loss 
rates and interactions energies in the two modes of the 
junction. Finally, we show that for strongly asymmetric 
losses one can protect the coherence of the MSGSs by 
suitably tuning the interaction energies, even after many 
loss events and in the absence of A^-measurement. 

Quench dynamics in a BJJ in the presence of atom 
losses. We consider a two-mode BJJ in the quantum 
regime. Initially, iVo atoms are in the ground state in the 
regime where tunneling dominates interactions, described 
by a spin coherent state 
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with 6 = tt/2 and = 0. Here is the joint eigen- 

state of the number operators fii in the mode i = 1,2 
(Fock state). The dynamics following a sudden quench 
of the tunnel energy to zero is given by the two-mode 
Bose-Hubbard Hamiltonian llSl 



ynj(n, - 1) 
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where Ei and Ui are the internal energy (or, for exter- 
nal BJJs, the lowest energy level of the well i) and the 
interaction energy between two atoms in the mode i, re- 
spectively, and Ui2 is the inter- mode interaction energy 
(for external BJJs U12 = 0). Setting h2 ~ No ~ fii, 
the interactions in Eq.(I2]) sum up to a non-linear term 
XTT-i, with X = (Ui + U2 - 2C/i2)/2. The time-evolved 
state at time tq = Tr/\xQ\ is a superposition of CSs, 

IV-^^Hi^)) = e-"^^«|V(0)) = ELUl^o;f,0fc), with 
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- 4>k = 27r/g and |cfe| = 
In the presence of two-body losses the markovian mas- 
ter equation for the density matrix p(t) of the condensed 
atoms reads [l9l \2(x 
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where 7^ and 712 are the loss rates of two atoms in the 
same mode i and of one atom in each mode, respectively, 
and {•, •} denotes the anticommutator. Since Eq.(j31) does 
not couple subspaces with different total atom numbers 
N = iii+h2 and N = Nq initially, at all times p{t) has the 
block structure p{t) ~ X^tvLo where WN(t) 
is the probability to have N atoms at time t and pN{t) is 
the corresponding conditional state. This decomposition 
is naturally^ accouiited for by quantum trajectories t i~>- 
li^jit)) = li^ji^)) IV'i(O) being the unnormahzed 
wave function when j loss events occur at times < si < 
52 < • • • < < i in the channels mi, • • • 
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For two-body losses one has three loss channels m = 
1. 2. and 12, the corresponding jump operators are 
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and M12 = 0.1(22. The dynam- 



ics between loss events is given by the effective Hamil- 
tonian = Hq — iD; the damping operator D = 
(1^1=1.2 ^ 1) + 7i2"i^2)/2 describes the gain of 
information on the system resulting from the knowledge 
that no loss occurred. The conditional state pN{t) after 
a detection of = A^o ^ 2j atoms at time t is obtained 
by averaging over the j jump times Sk and channels nik, 



PN{t) = WNit)pNit) = X] / ^-"^ 
mi,--- ,rnj 0<si < - - -<Sj : 
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where Pmu-,mjisi, . . . , Sj] j) = 7m, . ■ ■^ mAl'^ iHW is 
the joint distribution of the Sk, nik, and j [22|. By fur- 
ther summing over the number of jumps j one gets the 
total density matrix p{t) = X]jv=o '^N{t)pN{t) which is a 
solution of Eq.®. In order to understand the effect of 
losses we analyze separately each iV-atom sectors. 

Density matrix in the suhspace with Nq atoms. When 
no loss occurs in the time interval [0,t], from Eq.([4]) we 
obtain the unnormalized conditional state ^"^'^^ (t) = 
\ipo{t)){ipo{t)\ in the Fock basis. 
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where djVo("-i) = 5(7i+72-7i2)('t-i-"i)^ up to an irrel- 
evant constant, 2ni = [71— 72+^^0(272 — 712)]/ (71+72 - 
712), and p^°\t) = |?/>(°)(t))(?/>(°)(t)| is the lossless density 
matrix. For A^o ^ 1 the matrix elements of /5'"-'(i) have 
Gaussian moduli peaked at {ni, n[) = {No/2, No/2) with 
a width ~ VNq. For symmetric loss rates 71 = 72 and 
712 = 0, this peak coincides with the center of the Gaus- 
sian damping factor in Eq.®. Thus the MSGS formed 
at time tq is affected by damping when 71 > xi/^o- In 
contrast, for 72 = 712 = the damping factor is centered 
at (ni,ni) = (1/2, 1/2) and its effect on the MSGS sets 
in at the much smaller rate 71 « xq/^q- 

Density matrix in the suhspace with (Nq — 2) atoms. 
We focus now on trajectories t 1—^ |^i(0) having one 
loss event in channel m at the random time s G [0,t]. 
The instantaneous jump transforms a GS into a GS, 
Mm\No;e,(l)) cx |iVo - 2;6',(/«). This GS is rotated on 
the Bloch sphere by the evolution under the nonlinear 
effective Hamiltonian Hcs due to changes in energy and 
damping after the jump, yielding 



(7) 



where the random angles 9m{s) = 2arctan(e ^^'") and 
(^mis) = 2sXm depend on the random loss time s, the 
interaction energies xi = Ui - U12, X2 = -(C^2 - U12), 
X12 = {Ui — U2)/2, and the loss rate differences 61 = 
271 - 7i2, ^2 = -(272 - 712), and 612 = 7i - 72- Hence, 
apart from reducing N and producing damping, atom 
losses are identical to external 9 and 0-noises rotating 
the state around the z-axis by a complex angle 0m + 
ilntan(^) j2^. These noises have fluctuations S9m — 
\Sr,i\min{t,Ssm} and 6(j)m = 2|xm| min{t, 5s„i}, where 
6sm is the loss-time fluctuation. 

Tuning the energies at weak losses. Let us concen- 
trate on the MSGS formation time t = tq and take 
712 = 0. We first consider weak losses ^ Qx/Nq 
and symmetric energies Ui ~ U2. In this regime the 
0-noise is negligible since 69^ is much smaller than the 
quantum fluctuations ^1/ a/]Vo of a GS. In contrast, one 
has large ^-fluctuations S(j)„i ~ 2'!T/q, equal to the inter- 
component phase separation of the MSGS (we use here 
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FIG. 1. (Color online) Modulus of the density matrix ([8]) in 
the subspace with No — 2 atoms in the Fock basis at time 
t2, obtained by an exact diagonalization of Eq.([3]). Upper 
panels; symmetric losses rates (71 = 72 = x/(2007r)); lower 
panels: asymmetric losses (71 — 4x/(3007r), 72 = 0). Left 
column: symmetric energies (xi ~ — X2 = x)i right column: 
asymmetric energies (xi = 2x, X2 = 0). Other parameters: 
712 =Ei = 0, No = 100. 



Xi = — X2 = X Ssm > tq)- Such a strong (/)-noise has 
a relatively small effect on the coherences, albeit it trans- 
forms the corresponding statistical mixture of CSs into a 
mixture of Fock states (phase relaxation) [l^ . Actually, 
if 72 = the off-diagonal elements of the conditional den- 
sity matrix {ni,n2\pNo-2{tq)\n'i,n'r^ are weakly affected 
by losses for n'^ 7^ ni modulo 2, as shown in the ap- 
pendix and in Fig. [1] (lower left panel). The global shift 
to values ni, n\ < Nq/2 is due to the Gaussian damping 
caused by H^s in Eq.([7]). For 71 =72, instead, PNo-2{t2) 
is almost diagonal in the Fock basis due to a cancellation 
(occurring only for q = 2) when summing the contribu- 
tions of the two channels to = 1,2 (upper left panel in 
Fig[T]). We now consider asymmetric energies Ui 7^ U2- 
In order to keep tq = ir/lxql constant, we vary xij X2 
while fixing 2x = Xi ~ X2- Interestingly, it is possible 
to protect against 0- noise one channel, say to = 1, by 
choosing Xi = ^-nd X2 = ~2x, to the expense of en- 
larging noise in the other channel. Then 5(j)i ~ and 
S<f>2 = 4:71 /q. If only the first channel loses atoms, the 
conditional state after a single loss event is thus close 
to a MSGS with A^o ~ 2 atoms, apart from the damp- 
ing described in Eq.([n|). For symmetric losses, /5jVo-2(ig) 
has also large off-diagonal elements in the Fock basis due 
to probability one half to loose atoms in the protected 
channel (Fig. [1] right panels). 

The (/)-noise decreases at intermediate loss rates since 
the loss-time fluctuations Ssm decrease when increasing 



7,„. Indeed, we find Ssi w (27,; -^712) '^Nq ^ (i = 1, 2) for 
A^o 3> 1. Physically, at increasing j^n the loss has more 
chance to occur at small times, while for small 7™, s is 
equally distributed in [0,i]. 

Subspaces with Nq — 2j atoms. The wavefunction 
li^jit)) after j > 1 jumps is still given by Eq.Q upon 
replacing A'o — 2 by Aq — 2j and the angles of the GS 
by and ^ ELi^^M 



with tan(^^^ 



- ) . Thus the aforementioned effects per- 
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sist; for weak symmetric losses, though, keeping coher- 
ence by switching off phase noise in one channel is harder 
since the probability that all jumps occur in that chan- 
nel decreases exponentially with j . This means that for 
many loss events our proposal for protecting MSGSs is 
only efficient for strongly asymmetric rates. 

Total quantum correlations and sub-shot noise atom 
interferometry. We now show that for completely asym- 
metric losses not only the conditional states but also 
the quantum correlations in the full density matrix p{tq) 
can be protected by tuning the interactions. We mea- 
sure the amount of quantum correlations with the quan- 
tum Fisher information F{p), related to the best achiev- 
able phase precision in a Mach-Zehnder interferometer 
using p as input state by (A93)bcst — '^1 \/F{p) [2^ . 
Hence F{p) > (N) implies phase accuracy beyond the 
shot noise limit (Aiy9)sN = (N)^^^^ [25]. In our case 
Ftot(t) = F[p{t)] = ENlo^N{t)FN{t) where F^(t) is 
the Fisher information of the conditional state pAr(i) in 
the subspace with A^ atoms [2^. We optimize Ftot(t) 
over all directions of the interferometer. In a lossless 
BJJ, the two-component superposition has the highest 
possible value Ftot(i2) = A"q, nearly twice larger than 
that of MSGSs with q>2 components (13, 

Fig. [2] shows the Fisher information F]\f{t2) in the 
N — Nq ~ 2 subspace, corresponding to the state con- 
ditioned to a single loss event. For symmetric energies, 
the low values of Fn at small symmetric loss rates are 
direct consequences of the cancellation among channels. 
For 72 = a much larger value is found, showing that 
the aforementioned non-vanishing inter-component co- 
herences carry useful quantum correlations. At interme- 
diate rates, i^7v(^2) increases and reaches a maximum as 
a result of the reduced phase noise when 71 = 72 [s^- 
This reduction is clearly seen on the Husimi distribu- 
tions Qn{0,4>) = ^{N]e,(j)\pN{t2)\N]B,(l)) which display 
a flat profile for small losses and two emerging peaks for 
larger losses (insets of Fig. For asymmetric losses 
72 = 0, by choosing xi = we have FN{t2) — > A^^ in 
the small loss limit and the Husimi function has two 
peaks at the GS phases (p ~ ±7i'/2, in agreement with 
our prediction that the conditional state converges to 
a two-component superposition. At intermediate rates 
F/v(t2) decreases with 71 because the Gaussian damping 
compensates phase noise reduction. At large losses this 
damping transforms PNo-2{t2) into a superposition of 
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FIG. 2. (Color online) Fisher information in the subspace 
with A'^o — 2 atoms as a function of the loss rate 71 (in semi- 
log scale) at time t2 — T/A, with T = 2n/x- Solid black line: 
71 = 72, Xi = -X2 = X\ dotted red line: 72 = 0, xi = -X2 = 
X; dashed blue line: 72 = Xi = Oi X2 = ~2x; dot-dashed 
brown line: 71 = 72, Xi = 0, xa = — 2x- The optimization 
over the interferometer directions is done independently in the 
(A'^o — 2)-subspace. Insets: Husimi functions for the values of 
7i and Xi indicated on the curves by black circles. Other 
parameters as in Fig. [1] 
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FIG. 3. (Color online) Total quantum Fisher information 
Ftot{i) vs time t (in units of T). From top to bottom: 
71T = 8/300, 72 = Xi = 0, X2 = -2x (blue dashed line) 
7ir = 8/300, 72 = 0, xi = -X2 = X (red dotted line 
7iT = I2T = 1/100, xi = 0, X2 = -2x (brown dot-dashed 
line); 71T = 72r = 1/100, xi = -X2 = X (black solid line). 
Other parameters as in Fig[T] Dashed horizontal line: F%ot for 
the highest squeezed state in the lossless case. Inset: average 
number of atoms vs. time for the same parameters. 



Fock states with ni = or 1 atoms in the first mode [17| . 
For symmetric losses and even A'o , PNo-2{t2) is trans- 
formed instead into the Fock state | t) ti^-ving a larger 
Fisher information ^ 7V^/2. 

Our main result is presented in Fig. [3J which displays 
the total Fisher information i^tot {i) obtained from an ex- 
act diagonalization of Eq. ([3]) for various rates and ener- 
gies, keeping the same value of x ^^nd of the mean num- 
ber of atoms at time ^2 (see the inset). We find that 
for 72 = 0, tuning the energies has a strong effect on 



— 1/2 

Ftot{t) at times xt ^ ^0 > yielding to larger Fisher 
informations than for squeezed states. Especially for the 
two-component superposition, quantum correlations are 
then preserved even after the loss of 20% of atoms. 

Conclusions. We have shown that the interplay of 
atom losses and interactions in BJJs leads to different 
decoherence scenarios even after a single loss event. In 
particular, for strongly asymmetric two-body losses one 
can protect superpositions by tuning the interaction en- 
ergies, leading to useful states for high-precision atom in- 
terferometry applications. Such asymmetric losses occur 
in internal BJJs with ®^Rb atoms in the hyperfine states 
iFjUip) = |2,— 1) and |1,1), subject respectively to fast 
two-body and slow one-body losses Finally, we note 
that (/)-noise associated to the loss times Sk could be used 
to obtain an indirect measurement of Sk from the inter- 
ferometric estimation of the phase rotation of MSCSs, 
yielding to an observation of quantum jumps like in cav- 
ity QED [3. 

Appendix. Using Eq.([4]) we find the distribu- 
tion pi\s; 1) = 7„e-'^«'"[cosh(sJ„)]^«-2||e-**^"«|Afo - 



7m e 

2;ern{s),0)fNo{No - l)/4 with G, = (27, + 
27i - 2712 and G12 = (71+72 + 712)^^0 - 271 - 
Averaging over trajectories as in Eq.([5]) yields 

{ni,n2\pNo-2{t)\n[,n'2) oc 

^Cm{t;nl,n[) {ni , n2\p^^°[°2''^ (t) \ n[ 



^7r, 



712)^0 - 
272-712- 



(8) 



where p^^°_'°2''^(i) 



-itH, 



«|iVo-2;f,0)(7Vo"2;f,0|e'*^off 



is the conditional state having no loss in [0, t] for an initial 
CS with = iVo - 2 atoms, and 



Cm{t;n,n') = 



1 



-t[G„,+5,„(n+n'-Wo+2)+2iXm(n-"')] 



G,„ + 5.,n {n + n'-No + 2) + 2ixm{n - 71') ' 



Eq.® shows that the density matrix conditioned to a 
single loss event is a MSGS with A^o — 2 atoms modu- 
lated by the envelope ^ 7mCm(t; n, n') and the damp- 
ing factor of Eq.®. For weak losses j„i ^ Qx/^o and 



j: C/i C/2, |C™(i,;n,n')| ^ 3^|sinc( 



"h i vanishes at 
q n 

the values n' = n + pq {p = ±1,±2,---) correspond- 
ing to the nonzero off-diagonal elements of the statisti- 
cal mixture of CSs J2k \ck\'^\N] ^,(t)k){N] ^,<i)k\, and de- 
cays like |n' — n|~^ for the other off-diagonal elements 
encoding the coherences among the CSs [29|. Hence the 
coherences are non vanishing, except for 71 = 72 and 
5 = 2, because then the contributions of the two channels 
cancel each other, T™^™(^2; 't-, 7^') — for n ^ n' . 
Taking now Xi = and X2 = — 2x, one infers from 

Z]m7mGm(i2;?T-,?T-') ~ ^(71 +72'5„„') that PNa-l{t2) OC 

f}'p{°-2^\t2) when 72 = 0. Therefore, the conditional 
state is a MSGS slightly modified by Gaussian damping. 
The above calculation can be generalized to trajectories 
with j > 1 jumps between and t. For 'jmt ^ 1 and 
j ^ No, PNo-2j{t) is still given by Eq.([5]) upon replacing 
J2lmCm{t;ni,n[) by Em 7mC'm(i; "i, "i)]^ • 
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matrix (ni,n2\p'''^\tq)\n{,n'2) in the Fock basis is the sum 
of two matrices, the first one having zero matrix elements 
for n'l 7^ ni modulo q and corresponding to the statistical 
mixture of CSs, and the second one having zero matrix 
elements for n'l = ni modulo q and corresponding to the 
coherences I^fe^j,, Cfc4, |A; |, 0fe)(A; |, (^^z |. The latter 
do almost not contribute to the Husimi function but are 
responsible for the high values of the Fisher information. 
Note, however, that the total Fisher information Ftot{tq) 
decreases by increasing 7™ due to the rapid decay of the 
probability WNQ-2{tq) of losing only two atoms. 



